On the glueball spectrum of walking backgrounds from wrapped-D5 gravity duals 
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We compute the mass spectrum of glueball excitations of a special class of strongly-coupled field 
theories via their type-IIB supergravity dual. We focus on two subclasses of backgrounds, which 
have different UV-asymptotics, but both of which exhibit walking behavior, in the weak sense that 
the gauge coupling of the dual field theory exhibits a quasi-constant behavior at strong coupling over 
a range of energies, before diverging in the deep IR. We improve on earlier calculations, by making 
use of the fully rigorous treatment of the 5-dimensional consistent truncation, including the rigorous 
form of the boundary conditions. In both cases there is a parametrically light scalar glueball. In 
the first case, this is a physical state, while in the second case this result is unphysical, since the 
presence of higher-order operators in the dual field theory makes the whole (physical) spectrum 
depend explicitly on a (unphysical) UV-cutoff scale. 
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I. INTRODUCTION 

The study of strongly-coupled models of electroweak symmetry breaking (technicolor) is notoriously difficult, be- 
cause of the strongly-coupled dynamics itself [HQ- In particular, it has been the matter of a long debate in the litera- 
ture whether in a special class of such models, usually referred to as walking Q, a parametrically light pseudo-dilaton 
exists in the spectrum the phenomenology of which would resemble that of an elementar y H iggs particle P-Q. 
With the discovery at the LHC of a scalar particle with mass of approximately 126 GeV [i^. |ld|. this fundamental 
field-theory question has become of utmost importance also from the phenomenology perspective. If one could prove 
that no such light dilaton exists in any reasonable strongly-coupled theory, this theoretical fact, in combination with 
the experimental results, might suggest that electroweak symmetry breaking emerges in nature from a weakly-coupled 
sector. Vice-versa, purely on the basis of the LHC data, a number of phenomenological analyses [ll[ show that while 
the 126 GeV scalar is compatible with the elementary Higgs of the Standard Model, it could also be a composite 
dilaton emerging from stro ngly -coupled dynamics. 

Gauge-gravity dualities [12| Il3j offer an unprecedented opportunity to investigate the strong dynamics (a vast 
amount of lattice studies also appeared in recent years [HI): because they allow to rewrite untreatable strongly- 
coupled field theories in terms of extra-dimensional theories of quantum gravity in a more accessible limit in which 
calculations are possible. In particular, a completely algorithmic procedure exists for computing the spectrum of 
glueballs of strongly-coupled field theories, provided one knows not only the dual (10-dimensional) background, but 
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also the 5-dimensional sigma-model action obtained by consistent truncation |15l |lq |. This algorithmic procedure 
makes crucial use of the diffeomorphism invariance of the 5-dimensional sigma-model action, which allows to write 
the linearized equations for the fluctuations of the background in terms of physical (gauge-invariant) fields, with 
well-defined boundary conditions. Hence, it is nowadays possible to ask, within the framework of gauge-gravity duals, 
whether there are strongly-coupled, confining, field theories which admit a parametrically light state in the spectrum, 
to be identified as a dilaton. 

The most challenging part of this program is to identify suitable candidate models. It is known that within the 
effective field theory (EFT) framework (supplemented by naive dimensional analysis), it is very difficult to write 
models in which parametrically light scalars exist, because of UV-sensitivity at the loop-level. In principle, there is 
no reason why such a scalar could not be light, but in practice this requires fine-tuning to allow for cancellations 
between large independent contributions (for a radically different and unconventional perspective on this point, see 
for instance [l7|). This is at the root of the hierarchy problem in the context of models of electro- weak symmetry 
breaking, but the same types of EFT arguments (and open problems) appear also in other contexts, for example 
in inflationary cosmology. The mass of the light dilaton is not affected by such arguments in a fundamental (UV- 
complete) theory, in which the low-energy measurable quantities are not UV-sensitive. But then, one has to decide 
what kind of fundamental theory to study. The physics lesson one learns from the EFT arguments is that the generic 
field theory is unlikely to produce a light dilaton, the existence of which would emerge only in very special classes of 
dynamical models. For this reason the number of candidate models (in the fully rigorous context of top-down approach 
to gauge-gravity dualities) is very limited, in spite of many studies [l8l - l2lT ] . To the best of our knowledge, at present 
only one class of such models has been identified and studied in enough details, and [22| reported on the presence of 
a parametrically light scalar state in the spectrum. It is not known yet whether this state is a dilaton, nevertheless 
this result is of crucial importance: it can be thought of as a proof of principle, stating that there exist non-trivial, 
confining strongly-coupled gauge theories in which a parametrically light scalar emerges from the dynamics. 

A recent paper [23| could not find any light composite states. The reader should be alerted of the fact that the 
backgrounds considered in [23[ are those in [2l| (we will call these backgrounds class 2 in the main body of the paper) , 
which are different from those in [22j (which we will refer to as class 1). In this paper, we update the analysis of [22| |. 
refining the treatment of the boundary conditions and making it rigorous, along the lines of (l6| . We confirm the 
main results in 22] , showing that backgrounds of class 1 yield a parametrically light composite scalar in the spectrum. 
We also perform the calculation in [23j . analyzing some special backgrounds of class 2, by using the rigorous boundary 
conditions, and without any of the approximations used in [23j . We show not only that a light scalar is actually 
present in the spectrum also in class 2 (at least in some region of parameter space), but also that calculations in 
models of class 2 are difficult to interpret, since the whole spectrum of fluctuations is UV-cutoff dependent. In the 
limit in which the UV-cutoff is removed, the theory has a continuous spectrum. And on the other hand, keeping a 
finite cutoff means that solutions of this class admit an interpretation in terms of a four-dimensional field theory that 
is at best an effective field theory, for which one expects radiative corrections to recreate the fine-tuning problems. 
None of these problems arise in the first class of solutions, studied first in [22[ and for which we refine the analysis in 
this paper. 

The conclusion is that, after imposing the more rigorous (and easier to implement) boundary conditions of (l6j, w e 
confirm the existence of a light composite scalar, in important regions of parameter space, as originally stated in [22j j . 
We also discuss several other open problems and possibilities, which we will address elsewhere. 

II. GENERALITIES 

All the backgrounds of interest in this paper can be thought of as being generated by the strong-coupling limit of 
the background generated in type-IIB by a stack of N c D5-branes wrapping a special 2-cycle inside the base of the 
conifold, and hence are generalizations of the linear-dilaton solution in 24] . The five-dimensional description of the 
system is a subtruncation of the Papadopoulos-Tseytlin 12511 truncation, itself a subtruncation of the most general 
consistent truncation of Type-IIB related to the conifold [26| . 

The coordinates are labelled by the Minkowski the radial direction r (the only one that the ansatz will depend 
explicitly upon) and five angles with ranges < 9, 8 < n , < (j), <f> < 2n , < < Air. The PT ansatz [25[ (setting 
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a'g s = 1) is given by the following (Einstein frame) Type-IIB background functions: 

ds 2 = e 2 P- x ds? + {e x+ S + a 2 e x -s){ei + e 2 ) + e x ~~s (e§ + e 2 - 2a( ei e 3 + e 2 e 4 )) + e~ 6p - x el (1) 

dsl = dr 2 + e 2A dx 2 13 , (2) 

F 3 = N [e 5 A (e 4 A e 3 + e 2 A e\ - 6(e 4 A e 4 - e 3 A e 2 )) + dr A (<9 r 6(e 4 A e 2 + e 3 A ei))] , (3) 

H 3 = -h 2 e 5 A (e 4 A e 2 + e 3 A e 4 ) + dr A d r h Y (e 4 A e 3 + e 2 A e 4 ) - (4) 

9 r /i 2 (e 4 A ei - e 3 A e 2 ) + <9 r x(-e 4 A e 3 + e 2 A ei) , (5) 

P 5 = h + *^5, A = -/Cei A e 2 A e 3 A e 4 A e 5 , (6) 



where 



ei = — sin^ d0 , (7) 

e 2 = d0, (8) 

e 3 = cos "0 sin 6* g?</> — sini/; d9 , (9) 

e 4 = sin^sin^ c?0 + cos^ d# , (10) 

65 = dt/j + cos 8 d<fi + cos 9 d<fi , (11) 



together with the constraints 



/C = M + 2N(h 1 + bh 2 ) , (12) 

(e 2 9 + 2a 2 + e^a 4 - e^)dhi + 2a(l - e' 2 ^ + a 2 e~ 2 ~ 9 )dh 2 
X ~ e 2 9 + (l-a 2 ) 2 e- 2 9 + 2a 2 ' ( ' 

All other fields of 10-dimensional Type IIB supergravity vanish on the background. The constants M and N control 
the fluxes of the F$ and F 3 form, respectively 

We are interested only in solutions to the wrapped-D5 system, hence we set M = = hi = h 2 = K. = %, which 
leaves only the dilaton (f>, the metric and the F 3 form as non-trivial background functions. In the following we pick 
N c = AN = 1. The background functions (g,p, x, <f>, a, b) depend only on the radial direction r. The solutions we are 
interested in can be found by first making a change of variable dr = 2e~ 4p dp, and then putting [l9T | 

25= P 2 ~Q 2 c - 24 ^ ^(P 2 -Q 2 )P' 3 sinh(2p) 
(Pcoth(2p) - Q) 2 ' 6 131072 
c sx_ ^°(P 2 -g 2 ) 3 sinh 2 (2p) e4( ^_ sinh 2 2p 



8192P' ' 8(P 2 -Q 2 )P' (u) 

P 2p 



sinh(2p) (P coth(2p) -QY sinh(2p) 
_ e^°(P 2 -Q 2 )smh 2 (2p) 



e 



256 
where 

Q = 2/9C0th(2p) - 1 . (15) 

Here, we set the end-of-space in the IR at po = 0, and in order to avoid a nasty singularity in the IR we fine-tuned 
Qo = — N c (the general form of the solution can be found in (27|, for example). 

All the BPS equations of motion are assured to be satisfied provided P solves the (master) differential equation [13] 

P " + P '(^TW + ^ ± f-"»"'( 2 p))=o. W 

The most general solution to the BPS equations is controlled by 5 integration constants. We fixed two of them 
already, and we will comment in due time about <f> , which appears in the 10-dimensional dilaton. The last two 
integration constants parameterize the general P. 

The solution P — 2N c p gives the background of [Hj]. Notice how this solution depends only on N c , but there are 
no integration constants. There are two main classes of solutions, with different types of UV-asymptotic behavior. 
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• Class 1: backgrounds where P is approximately linear with p at asymptotically large p. The solution in [24[ 
belongs to this class, together with the class of walking backgrounds studied in [22| . The field-theory interpre- 
tation of these solutions is explained for example in [19J , and is connected with the infinite Higgsing of the dual 
quiver gauge group. The dilaton is linear (in the UV). This class is characterized by one integration constant 
only, namely the value of -P(O) = Pq. 

• Class 2: backgrounds where P ~ e 4p / 3 diverges exponentially in the far-UV. The solutions discussed in [2ll an d 
studied by [23| belong to this class. As explained for example in (by drawing heavily on results from j28|). 
the dual field theory of these solution contains a dimension-8 operator, deforming the theory and making the 
UV badly behaved. The dilaton is approximately constant in the far-UV. These solutions are characterized by 
two integration constants. 

The Papadopoulos-Tseytlin [25| ansatz can be rewritten by promoting all the functions appearing in the background 
to eight five-dimensional scalar fields $ a , with sigma- model action defined by [l5[ 



S = J d b x^—g (Ir - lG ab d M <S> a d M $ b - U($) ) . ! 17 ! 



with kinetic terms 



G ab d M <^ a d N ^ b = \d M gd N g + d M xd N x + 6d M pd N p + ^d M <t>d N <£ (18) 
+ ^e- 2 ~ 9 d M ad N a + ~N 2 e*- 2x d M bd N b 

— <& — 2x 

(l + 2e~ 2S a 2 )d M hid N h 1 



e 2 ~9 + 2a 2 + e- 2 ~9(l-a 2 ) 2 
1 
2 



+ -(e 2 ~ 9 + 2a 2 + e~ 2 ~ g (l + a 2 ) 2 )dMh 2 d N h 2 + 2a(e- 2S (a 2 + 1) + l)d M hid N h 



and where the potential is 



V = -^e 2p - 2x (e 9 + {l + a 2 )e- 9 ) (19) 
8 

+ ±a 2 e~ 29+Sp 



l iY 2 e *-2x+8p [ e 2g + £ -2g {a 2 _ 2ab + y2 + 2(fl _ ft) 2 

8 



1 

r 



-4>-2x+8p^2 



2 



+ le 8p - 4x (M + 2N{h x + bh 2 )) 2 . 
8 

Again, in our case h± — h 2 — M = means that the last two lines of both the kinetic terms and the potential drop. 
The full five-dimensional PT system yields a vast space of possible background configurations, in principle controlled 
by 16 integration constants. This is a consistent truncation of Type-IIB, in the sense that given any such background, 
one can construct the full 10-dimensional background in a purely algebraic way, and what results is a solution of 
the classical equations of motion of Type-IIB. And the wrapped-D5 system (with six scalars) is itself a consistent 
truncation of PT. 



A. Fluctuations 



In order to compute the spectrum of scalar glueballs of the dual theory, one can study the fluctuations of the 
background in the five-dimensional sigma-model. This will not yield the full spectrum, but a subsector of it, which is 
going to contain the most important states at low energies. In order to do so, one has to fluctuate all the scalars $ a 
remaining after truncating (six in our case), but also the five-dimensional metric. We use the conventional ordering 
$ a = {g,p, x, </>, a, b}. In doing so, one has to keep into account the fact that S is diffeomorphism invariant, and 
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hence many of the fluctuations are unphysical (pure gauge) . The problem is overcome by making use of appropriate 
combinations of the original fluctuations, constructed in such a way as to preserve gauge invariance in a manifest 
way. These are called gauge invariant variables. Following [TBI ], it turns out that the number of such propagating 
gauge- invariant fluctuations a° is the same as the number of independent scalars in S (again, six for present purposes). 
The fluctuations a a satisfy the following equation of motion in the bulk [lf| 



Vl + dA'V r 



with boundary conditions 



e - 2 V 



V a ]c - K a , 



bed 



4(<S>' a V c + V a <5>' c ) 16V<S>' a <f>' c 



(d-l)A' 



(d - I) 2 A' 2 



0, 



(20) 



5% + e 2A q-' 2 (v a - dA'<S>' a - A Q | C $ ,C 

2 



2$I 



V r a L 



2&°<f>> b 
l fc+ (d-l)A> 



e 2A q~ 2 



{d-l)A 



(d - I) A' 
(v a - dA 1 ®" 1 - A°,,$' c 



4y<f>' b 

(d- 1)4' 



(21) 



The details about the derivation of these results can be found in [l5|,[l6[. Here, we only summarize the meaning of the 
notation: D r is the covariant derivative with respect to the sigma- model metric G ab , d = 4, q 2 relates to the Minkowski 
4- momentum with signature {+,—,—,—}, Tt a bcd is the Riemann tensor associated with the sigma-model metric G abl 
which is also used to raise and lower the sigma-model indexes a — 1, ■ • • ,6, V a — dV/d<& a , primed quantities refer 
to derivatives in respect to r, and all the functions A, TZ, $ a , V and their derivatives are evaluated on the classical 
background. We refer the reader to [HI, [l6| for more details about the formalism. 

Notice how the boundary conditions depend on undetermined coefficients A a | fe the meaning of which is explained 



at length in [Iff. In the limit of A a , 



\b\ 



-e 2A q- 2 



±oo, which we adopt from now on, we obtain the boundary conditions 



2$'' 



(d - 1)A> 



(d - 1)A' 



Vb 









a" 


= a a 








r 



(22) 



The spectrum can be computed by solving the bulk equations and imposing the UV and IR boundary conditions at 
two arbitrary end-of-space points in the radial direction, that we denote r± <C T2- These serve as regulators. The 
calculation has to be repeated for larger values of r-i and smaller values of T\ approaching the physical end of space. 
If the background is sufficiently well-behaved, this procedure converges, and the final results do not depend on the 
regulators. 

After changing the radial coordinate from r to p, we can write Eq. (f20f as (for k = — A — 4p + log 2) 



S a b d 2 + S a b d p + T a b + S a b e 2k q 2 



a b = 0, 



(23) 



where 



S a b =2G\ c d p <S> c + 4 (d pP + d p A) 5 a b , 



n =d b G a cd d p ^ c d p ^ d - 4e" 8p 



f4(V a d p ^ c + d p <i> a V c ) 16Vd p ^> a d p ^ , „ 
3^4 + 9(d p A) 2 1 Gcb + 9bV 



We will see that this further rewriting is going to facilitate the analysis 
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III. WALKING BACKGROUNDS WITH LINEAR-DILATON ASYMPTOTICS 

We start from the walking backgrounds of class 1. The spectrum has already been studied in (22j. However, while 
making use of the gauge-invariant bulk equations, in that paper the boundary conditions have been chosen on the 
basis of regularity arguments, because the general form of the boundary conditions was not known at the time. We 
are now in the position to redo this calculation, imposing the rigorous boundary conditions. 

We proceed as follows. First of all, we choose a large class of solutions P of the master equation which have 
linear UV-asymptotics (a few examples are shown in Fig. [IJ. All of them differ only by the value of P(0) = Po, or, 
equivalently, by the value p„ of the coordinate at which the linear behavior appears, taking over from the approximately 
constant behavior of P in the deep IR. Then we introduce two cutoffs p\ and p2, and fix <p a = 0. We compute the 
spectrum for all the backgrounds. We repeat this for various choices of p2 and pi, until the results show no dependence 



L , , , 1 , , , 1 , , , 1 , , , 1 , , , U 0.0 

2 4 6 8 10 02468 10 12 

P Po/2 

FIG. 1: The left panel shows P for a few different backgrounds. The right panel shows the spectrum q 2 computed for different 
values of Po/2 ~ p*. As Pa — ¥ 0, the spectrum agrees with the one for the solution in [24|, as computed in [15l ]. For large Po, 
there is a light scalar. These numerical results have been obtained with p2 = 300, and we checked that different (but large) 
values of this UV cutoff do not appreciably change the results. 
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Po/2 Po/2 

FIG. 2: The left panel shows - q 2 versus Po/2 ~ p*. The right panel shows a log- log plot of the mass squared of the 

lightest scalar versus Po/2 ~ p„, with the blue line corresponding to q 2 — 2/Pq. 



on the cutoffs. The result is shown in Fig. [T] For completeness, we also show some more details about the spectrum 
in Fig. H 

Let us discuss the results. First of all, the spectrum is discrete only up to q = 1. For larger values of q the system 
of linear equations breaks down: while there are discrete solutions for finite p2, taking p2 — > oo yields a continuous 
spectrum. In order to understand where this comes from, let us expand in the far-UV: 



P = 2p 

Neglecting exponentially suppressed terms, this yields 



0(e 



-4p\ 



s a h 
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o( P - 



T1 



( 


-4 





4 


-2 





° \ 





-6 


-1 


l 

2 










2 


-6 


-3 


1 

2 










-4 


-12 


2 


-3 






















-4 


4 


\ 














4 


-4/ 



So that all the non-zero entries of these matrices are 0(1) in the UV. But we also find that 



(25) 



0i 



(26) 



e 2fc - 1 + O (e- ip ) 



(27) 



After going to a diagonal basis, one finds that the fluctuations behave (up to powers of p) as 



with C = 1 for half of them and C — 9 for the other half. This is the origin of the problem: we would like to 
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interpret the (real) solutions of the bulk equations in terms of operators of the dual theory and their sources, but 
this cannot be done for large-g 2 because there is no sense in which half of the asymptotic solutions are dominant 
and half subdominant, rather all the solutions become oscillatory in the far UV at large-g. This is a well known 
behavior that appears in many quantum mechanics systems: provided one keeps the UV cutoff as physical, there is a 
discrete spectrum (similar to what happens for a potential well of finite size in quantum mechanics) , but if the cutoff 
is removed, the spectrum becomes continuous. The result is similar, qualitatively, to the hydrogen atom in quantum 
mechanics: a discrete spectrum at low energies, with the eigenvalues becoming closer and closer to each other until a 
critical value beyond which the spectrum is continuous. We adopt the same interpretation: above q = 1 the states are 
unbounded. The existence of a theory in which confinement yields a spectrum of bound states with an upper bound 
on the mass is at least unconventional. But we do not need to worry about this: the low-energy behavior is consistent 
with expectations. 

Looking at the low-g 2 spectrum of excitations yields a very interesting result: there is a parametrically light state, 
the existence of which emerges only when p* is large. The precise numerical values of the masses are slightly different 
from those in [22| but qualitatively the same. This state does not exist for the classical linear dilaton solution [24| . 
and indeed our results agree (within numerical accuracy) with those of 1 1 51) in the p* — > limit. All of this means 
that we have obtained results which are in qualitative agreement with [22j |. 



IV. WALKING BACKGROUNDS WITH CONSTANT-DILATON ASYMPTOTICS 

The walking backgrounds of class 2 are very different. In the far-UV the exponential behavior of P ~ e 4/ V 3 is the 
consequence of the fact that in the dual field theory there is a dimension-8 operator. The subtle relation between 
this fact and the baryonic VEV, and the fact that it is possible to UV-complete the dual theory, by algebraically 
constructing new solutions of type IIB in which F$ and E>2 are non-trivial, is explained in [L9| . We will not further 
pursue this line, which requires a dedicated study. 

Our attention will focus on the UV-cutoff dependence of the glueball spectrum, and on the effect of the rigorous 
boundary conditions. To do so, we start from the UV-expansion of the relevant quantities appearing in the bulk 
equations for the fluctuations. In the UV, these solutions behave as 



P = 3c+e 4 "/ 3 + — I p 



4 

3cj 



13 

> + 16 



,-4p/3 



8c + p 



192c 2 + 



,-Sp/3 



01 



(28) 



The two integration constants c+ and c_ characterize completely the background. The former is related, on the field 
theory side of the duality, to the insertion of a dimension-8 operator, the latter to the appearance of a dimension-6 
VEV which dominates the IR physics (which is ultimately the reason why one might expect a light dilaton in the 
spectrum) . For illustration purposes we show in Fig. [3] a sample of possible solutions P obtained by varying c + and 
c_, or equivalently Pq and p». 

With these expansions, we find that 
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+ O (e- 4 "/ 3 ) , (29) 



in which all the non-zero entries are 0(1). However, the factor multiplying q 2 in Eq. (|23|) behaves as 



'2k 



= 2c + e 



4p/3 



o 



(30) 



growing without bound in the UV. This means that for any value of q 2 , the fluctuations become oscillatory in the 
far UV, so that the spectrum only consists of a continuum, and therefore, contrary to what happens in class 1, there 
is no range of q 2 in which a discrete spectrum can be obtained. The only way to get a discrete spectrum of bound 
states is by assuming that a UV-cutoff is physical. But then the results of the calculations will depend explicitly on 
this cutoff, and hence they cannot be interpreted as the spectrum of a UV-complete strongly-coupled confining field 
theory. As anticipated, a possible way to provide a UV-completion exists, and the study of the spectrum in this case 
will be presented elsewhere. 

For completeness, we computed the spectrum for three particular solutions in this class, by varying the value used 
for the cutoff. We report the results in Fig. 2) Notice how the whole spectrum depends strongly on pi when p2 is very 
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P 




1 u — 1 ' ' — 1 ' ' — 1 — 1 ' ' — 1 — 1 ' — 1 — 1 — 1 ' — 1 — 1 — u 10 u ' — 1 — 1 — 1 1 ' ' — 1 ' — 1 — 1 — 1 1 ' — 1 — 1 ' — 1 — 1 — u 

02468 10 02468 10 

P P 

FIG. 3: The function P, solution to the master equation, for various examples of solutions in class 2 (notice the logarithmic 
scale). In the left panel a subclass having all the same leading asymptotic behavior in the UV, but different values of Pq = P(0). 
In the right panel, a set of solutions for P having the same value of Pq, but different UV-asymptotics. 



large, to show explicitly that the discrete spectrum of bound states actually does not exists, but is an artifact of the 
UV-regulator. But notice also that if one interprets p2 as a physical scale, and keeps it fixed (and smallish), the result 
is a spectrum in which one state is anomalously light, at least is some region of the parameter space. Notice also that 
the numerical value of the light mass obtained in this way, for the examples (a) and (b) in Fig. [4] agrees reasonably 
well with the mass computed for the same value of Pq in Fig. [TJ but then this result looses meaning when taking p2 
large, because the state is overwelmed by the fact that the heavy states all become light. In the case of example (c), 
the situation is more subtle: for small values of pi is it not possible to really gauge whether there is a sense in which a 
light state is present. One important difference with [23j is in the fact that we used the rigorously derived boundary 
conditions from [l6j]. However, the existence or not of a light state in all of these backgrounds cannot be established 
here, but requires redoing the analysis for backgrounds in which the UV-asymptotic behavior is more tame fl9j | . 

There is also another significant difference between our analysis and the one in [23| . To understand this, we need 
to digress, and explain to the reader more in detail the meaning of the integration constants in P, in relation to the 
properties of the master equation. By looking at the master equation, and the relation between P and the functions 
appearing in the supergravity ansatz, one finds that obtaining a smooth solution requires imposing the requirement 
that P and P' be monotonically increasing, and satisfy always P > Q. By inspection, it turns out that actually one 
has to require P > 2N c p. Also by inspection, one finds that there are only three possible behaviors admissible for P: 
it can be approximately constant, linear with p, or exponential with P ~ e 4p / 3 . Together with the requirements of 
monotonicity, this means that any solution can be approximated by three sections over which P(p) is approximately 
constant at small-p, linear at intermediate-/), and exponential at large-p. Such is the behavior of example (a) studied 
and shown in Fig. 2J One or more of these behaviors may be absent, provided one always has P > 2N c p, such as for 
(b) and (c) in Fig. 2) The walking region is the one in which P is approximately constant, as explained in [21(. 

In [2l[ a special class of solution was identified, in which the linear behavior is absent. The general solution for P 
was shown to be written as a power expansion in a parameter c such that 

P = cPi + ip_i + ip_ 3 + • • • , (31) 

c c d 

and in Figs. 1, 2 and 3 of [2l| the goodness of such an approximation was gauged by comparing to the exact 
numerical solution of the master equation. The result is that for any physical quantity that does not depend on the 
10-dimensional dilaton <f> the agreement between the leading order approximation P ~ cP\ and the exact result is 
good, but the convergence of the series is slow for the dilaton, in which the leading-order approximation is visibly 
very different from the exact result, even at the qualitative level {<fi appears to be non-monotonic in the approximate 
expression). Hence, any calculation where one needs to use the dilaton requires using the numerical result, unless 
one dials to unreasonably large values the constant c, such that for all practical purposes the dilaton becomes exactly 
constant. In particular, the calculation of the spectrum of scalar excitations is one such delicate case, because the 
dilaton appears in a very non-trivial way both in the sigma-model metric and in the potential. 

The physical meaning of the parameter c in this expansion is closely related to the dimension-2 baryonic VEV. A 
complete discussion of this subtle issue can be found for example in (and relies on results from [H,[2^]). Large 
values of c are related (again, non-trivially) to taking small values of this VEV. At a more technical level, one finds 
that for a given Pi (p) there exists a minimum value of c such that the constraint P ^ 2N c p is satisfied for all p [2l[ : 



P-2 





FIG. 4: Spectrum obtained for examples of P in class 2. The upper-left panel shows the three choices of P we analyze, labelled 
as (a), (b) and (c). The other three panels show the spectrum of q 2 obtained for these solutions, as a function of p2. Example 
(a) is shown in the upper-right panel, example (b) in the lower-left and example (c) in the lower-right panel, notice that all 
the spectra degenerate to a continuum for large p2, and that in some cases an anomalously light state is present as long as p2 
is kept fixed and small. 



in this case, P becomes very close to the linear dilaton solution, at least for some intermediate region of p. The case 
of the solutions of class 1 is such that this minimum value is realized, in such a way that asymptotically P ~ 2N c p in 
the UV. On the contrary, large values of c are realized by solutions that are always very different from the solutions 
in [24] and (more relevant to this paper) those in [22} . 

All of this means that the regime in which the approximation P ~ cP\ can be used is the opposite region of 
parameter space in respect to the one in which the analysis in [22j identified a light composite scalar, and hence the 
analysis in (23l |. which adopts this approximation throughout the whole paper, cannot be compared to (22j. 



CONCLUSIONS 



In this paper we computed the spectrum of glueballs of the theories the gauge-gravity dual of which is represented 
by walking solutions to the wrapped- D5 system in Type-IIB. We considered two possible classes of solutions, which we 
name class 1 and class 2, and that differ by the UV-asymptotic behavior. We applied in full rigor the tools of consistent 
truncation and of the gauge invariant formalism developed in [llj , supplemented by the boundary conditions of [IF 

In class 1 we find qualitative agreement with the literature: in spite of the fact that the boundary conditions in 
had been chosen on the basis of regularity and superficially look very different from those constructed rigorously in 
and applied here, both analyses show that there is a parametrically light scalar state, the mass of which is suppressed 
by the length of the walking region. It is still to be proven that this is a dilaton, but because it is known that the 
walking region is related to the presence of a VEV for a dimension-6 operator in the dual theory, it is possible that 
this is actually the case. This question has to be investigated further. 

For class 2, we perform the analysis of the spectrum for a sample of possible solutions. First of all, because P 
diverges exponentially in the UV, and hence a higher-dimensional operator is present in the dual theory, spoiling its 
UV-completeness, it is not possible to remove the UV-cutoff from the calculation, and hence all the results depend 
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explicitly on such a cutoff, making their interpretation questionable. Secondly, even if one takes the UV cutoff to be 
fixed, a light state is present in the spectrum in the interesting range of parameter space where the comparison to 
backgrounds of class 1 can be made [22[ ■ 
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Appendix A: Goldberger-Wise mechanism and boundary conditions 

We consider here a simple exercise, as a commentary on the use of the boundary conditions in Eq. (|2"2"|) . We 
reconsider one of the examples in (l6| . namely the spectrum of the simplest possible realization of the Goldberger- 
Wise mechanism, already discussed and well-known in the literature ■ We do so in order to illustrate how dangerous 
it is to use naive choices of boundary conditions, in place of the rigorously inferred Eq. (j2"2)) . 

The model is defined by a sigma- model with only one scalar $, having canonical kinetic term, and superpotcntial 

W($) = -\-^ 2 , (Al) 

with the potential given by 

The background solution we are interested in is 

$ = $ie" Ar , (A3) 
A = r-i^e" 2 ^, (A4) 

and for simplicity we set r± = and keep it fixed. We will also keep r 2 ^> n fixed. 

Following 0,111, we take $i G 0(1), and assume that r<i ^> 1, and then solve the bulk equations for the fluctuations 
in the limit in which q 2 — 0. This yields 

a (r) = Cl e- Ar + c 2 e-( 4 - A > , (A5) 

with ci and c 2 arbitrary real integration constants. We then plug the result into the IR and UV boundary conditions 
and solve for c\/c2 and q 2 . The result is well known 0: 

^ 4A 2 (A-2)3>f 

9 ~ 3(-l + e2(A-2)r 2) - W 

Notice that in |16j this result was also tested numerically, without making use of any of the approximations mentioned 
above, obtaining excellent agreement for choices of parameters for which the approximations yielding this result hold. 

Some comment about the physical meaning of this result will be useful later. When A > 2, the background scalar 
represents a VEV in the dual, otherwise conformal, field theory. The presence of a UV cutoff is the only explicit 
source of breaking of scale invariance, and hence in this case the mass of this lightest state is parametrically small. 
In the case |A| <C 1, the scalar represents the insertion in the dual theory of a quasi-marginal operator. The explicit 
breaking of scale invariance is hence small, while the presence of the IR cutoff induces spontaneous breaking. Again, 
one finds a parametrically light scalar. The last case, when A < 2 is not small represents the generic case in which 
the effects of explicit symmetry breaking are larger than the effects of spontaneous breaking, and the resulting mass 
is not small. This is also the case in which the approximations we made break down (unless $i <C 1), as shown in pi 
for example. All of this is in perfect agreement with field-theory expectations. 
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Let us look more in detail at the practical effects of the UV boundary condition. We can write it explicitly (at 
leading order in the approximations we are making): 



3 



(2-2A)r 



(d + A) a 



(A7) 



Noticing that (d r + A)a = (—4 + 2A)c 2 e~( 4 ~ A ) r , one sees that in the limit in which ri becomes parametrically large, 
this boundary condition reduces to the requirement that the subdominant asymptotic solution is kept, while setting 
the dominant solution to zero. This would be the result also if we were to impose the simpler a{r2) = boundary 
conditions. The advantage of using Eq. (|2"21 is not clear in the case of one scalar and of a background which is known 
in closed form, but when many scalars are present, they mix, and the background is known only numerically, this is 
clearly a much more effective way of writing the boundary conditions. In particular, in the case of mixing between 
many scalars, it may be difficult to give a precise meaning to what it means to be dominant or subdominant for any of 
the possible asymptotic solutions, since in each scalar there will be more than two independent coefficients appearing 
at various orders (because of mixing, again). Also, this equivalence of choices is true only provided the r 2 — > +oo limit 
can be taken. If one wants to interpret T2 as a physical cutoff (kept fixed) the three choices yield different results, 
which converge only for large r^. We remind the reader that Eq. (|22[) descends directly from the requirement that 
the variational principle be well-defined for the full five-dimensional action, and hence this is the correct choice (up 
to the ambiguity in the coefficients mf discussed in and not relevant here). 

Let us now perform the exercise of computing again the mass of the lightest scalar state, but this time we will make 
use of an arbitrarily modified form of the IR boundary conditions. For simplicity, and considering that we have only 
one scalar for which a semi-analytical solution is known, we impose a(r2) = 0, following the indications discussed 
earlier, and restrict our attention to very large values of r 2 . The form of the IR boundary conditions we want to test 
is controlled by two arbitrary parameters a and /3, according to 



9 e (2-2A)r 

^A 2 d> (ad + /3 A) a 

6 q z 



(3 a 



(A8) 



With this choice, a = means Dirichlet boundary conditions, while /3 = means Neumann boundary conditions. 
After imposing the UV boundary conditions, and taking again T2 3> 1, we find 



2A 2 $ 2 / a 



2— - A + —(A - 2)cotanh ^ (A - 2)r 2 



(A9) 



For generic values of a and j3, we see that q 2 £ 0(1), hence failing to reproduce correctly the behavior expected for 
A > 2. By making use of the expansion cotanh(x) ~ 1 + 2e~ 2x we have (for A > 2) 



2(a-/3)A 3 $f 4a 



3/3 



A 2 (A- 2)$ 2 



-2(A-2)r 2 



(A10) 



Clearly, the first term is parametrically suppressed only for |A| -C 1, but not for A > 2. The second term is 
parametrically suppressed when A > 2, but in order to obtain the light state one must then cancel the first term, 
which requires a — /3, which is the choice of boundary conditions we adopted in this paper! Notice in particular that 
a = corresponds to Dirichlet boundary conditions, and in this case there cannot be light states unless A is very 
small. Even worse, taking /3 = (Neumann) yields to a singular behavior, which just means that there is no light 
state for any possible choice of the parameters. 

In conclusion, this simple little exercise shows that the boundary conditions we adopted in the main calculations of 
this and related papers, namely Eq. (|22|) . reproduce the physical results (as we already knew from @, EE])- Imposing 
Neumann or Dirichlet boundary conditions yields unphysical results. Ultimately, the reason for this is that the correct 
boundary conditions have been derived by imposing consistency with the variational problem. 
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